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Abstract 

The technique of differential intertwining operators (or Darboux transformation 
operators) is systematically applied to the one-dimensional Dirac equation. The fol- 
lowing aspects are investigated: factorization of a polynomial of Dirac Hamiltonians, 
quadratic supersymmetry, closed extension of transformation operators, chains of 
transformations, and finally particular cases of pseudoscalar and scalar potentials. 
The method is widely illustrated by numerous examples. 
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1 Introduction 



In recent times a growing interest to applications of supersymmetric quan- 
tum mechanics (SUSY QM) in different fields of theoretical and mathemati- 
cal physics is noticed. Thus, it has been applied to generate new families of 
exactly solvable potentials [l]-[4], to study ordinary symmetries of the set of 
Riccati equations [5], and also to analyse a new kind of symmetry of ordinary 
differential equations, called translational invariance with respect to Darboux 
transformations [6]. As it has been understood after the work by Andrianov 
et al. [7], SUSY QM (first introduced by Witten [8] as a toy model in Quan- 
tum Field Theory) is basically equivalent to a method of finding solutions 
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of a differential equation from the known solutions of another equation just 
by applying to them a differential operator. This approach was studied for 
the first time by Darboux [9] and is known nowadays as the method of Dar- 
boux transformations. Its most famous applications are related with nonlinear 
equations of mathematical physics [10] and inverse scattering problem [11]- 
[13]. Another fornmlation of the same method is due to Schrodinger and it is 
known as factorization method in quantum mechanics [14,15]. 

Though this method is widely used for Schrodinger-hke equations, its applica- 
tion to systems of equations, and in particular to the Dirac system, is studied 
much less. Probably the first paper where a differential- matrix intertwining op- 
erator for the Dirac equation is constructed is due to Anderson [16]. It deals 
with component equations for the transformation operator. For the particu- 
lar case of a pseudoscalar potential this author has found rather complicated 
equations defining the transformation operator. This might be the reason why 
his method has not received much popularity. In another paper [17] Stahlhofen 
applies a general theorem proved in [10] to study transparent pseudoscalar po- 
tentials for the Dirac equation, obtaining and analysing new kind of potentials 
with eigenvalues embedded in the continuous spectrum. Although there is not 
a new technique in this paper, it illustrates the idea that similar methods may 
find unexpected applications. We notice also the papers [18,19] (for a review 
see [10]) where ad hoc Darboux transformation operators are constructed for 
some systems of differential equations, and the paper [20] where the authors 
use differential transformation operators for finding particular solutions of the 
inverse problem for a regular Dirac operator. 

We would like to remark that the technique of integral intertwining operators 
is well-known in the method of inverse quantum scattering [21] (see also [22]), 
but difi'crcntial-matrix intertwiners have been found only quiet recently [20,23]. 

In this paper we give a consistent introduction to differential-matrix inter- 
twining operators for the one-dimensional stationary Dirac equation (Section 
2) and study their main properties. We shall show that every property of 
a differential intertwiner for the Schrodinger equation finds its counterpart 
in the case of the Dirac equation. Thus, in Section 3 we prove a theorem 
about factorization of a polynomial of Dirac operators by means of Darboux 
transformation operators, which serves as a basis for establishing a hidden 
quadratic supcrsymmctry underlying our method (Section 6). In Section 4 
we show that despite the fact that our transformation operators have non- 
trivial kernels, a one-to-one correspondence between the spaces of solutions 
of the initial and transformed equations is possible. Section 5 is devoted to 
analysing transformation operators as acting in a Hilbert space. In Section 7 
we study chains of transformations and establish relativistic analogs of Crum 
determinant formulas [24] . Sections 8 and 9 are devoted to studying particular 
cases of pseudoscalar and scalar potentials. Our method is widely illustrated 
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by numerous examples given in Section 10. Finally, in Section 11 we draw 
some conclusions and give an outlook of future works. 



2 Darboux transformation for the Dirac equation 

Consider the one-dimensional stationary Dirac equation 

ho^ = . (1) 

Here ho = jOx + Vq, 7 = ia2, 02 is a Pauli matrix, = djdx^ Vq is 2 x 2 real 
symmetric potential matrix and the function is a two-component vector (we 
shall call it spinor), "tp = (V'l, V'2)* (the superscript "t" meaning the transpose). 
We do not indicate the interval where the real variable x falls, because our 
constructions are independent of this interval. 

Suppose we know all solutions tjj to Eq. (1) and we want to solve a similar 
equation for another Hamiltonian hi defined by the potential Vi. The problem 
of finding eigenfunctions ip of the Hamiltonian hi , 

hi^ = E^, (2) 



can be reduced to the problem of finding a transformation operator (see e.g. 
[22]), also called intertwiner [5], defined by the equation 

Lho = hiL. (3) 



It follows from here that the function (/? = Lip is a solution (maybe trivial) 
of Eq. (2). Here and in the following we shall suppose that all operators act 
in the space of spinors with sufficiently smooth (i.e. infinitely differentiable) 
components. In this paper we restrict ourselves to differential transformation 
operators. From this point of view let us consider the simplest differential 
transformation operator 

L^Ad^ + B. (4) 



Here A and B are 2x2 matrices with x-dependent entries. Similar trans- 
formation operator for the Schrodinger equation leads (see e.g. [25]) to the 
transformation investigated for the first time by Darboux [9]), and therefore 
it is usually called "Darboux transformation operator". In this context it is 
natural to call any differential transformation operator Darboux transforma- 
tion operator [25]. 
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Using the intertwining relation (3), and assuming the hnear independence of 
derivative operators d^, k — 0, 1, . . ., we are led to the following system of 
equations for yl, B and Vi. 



A-f-jA = 0, (5) 
AVo - V,A + B^-^B-^A,^0, (6) 
AVo:, + BVo - 7^. -ViB^O. (7) 

The subscript x denotes the derivative with respect to x, e.g. A^ = dA/dx, 



and the matrix A^ is composed of the first derivatives of the elements of the 
matrix A. 

The Eq.(5) may be considered as a restriction imposed on the matrix A: two 
of its four elements may be chosen arbitrarily, and the other two should be 
found in terms of the previous ones. The Eq. (6) defines the potential difference 

D = {AVo - VoA + B^-^B- ^A,)A-' . (8) 

We are supposing that A is an invertible matrix for all values of x. Prom (7) we 
will find B. To integrate this equation we first change the dependent variable 
B in favor of a new variable fl. B — Afl. Thus, upon using (8) we obtain from 
(7) the equation for fl 

(Vo - 7ii)^ + nvo - VoQ + (7^1 - n-f)n = o. (9) 



This equation can be considered as the derivative of a relativistic analog of 
the Riccati equation that appears when similar calculations are carried out 
on the Schrodinger equation. Eq. (9) can be linearized and integrated by a 
similar substitution: ^2 = —UxU~^. After such a calculation we get 

[W-^(W + 7W,)], = 0, (10) 



where U is supposed to be everywhere invertible. Now, denoting by A a matrix 
integration constant, we get the equation for U: 

7W^ + yoW = WA. (11) 



The matrix A should be taken hermitian. 

We notice that (11) is very similar to the initial Dirac Eq. (1). The single 
difference is that ip is a. spinor and E is a. number, while U and A are matrices. 
Nevertheless, if solutions of (1) are known we can immediately solve (11). 
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Since the matrix A is hcrmitian it always can be reduced to a diagonal form. 
Therefore let us take A to be diagonal 



A 



^ Ai 
Ao 



(12) 



If we now choose U — (1*1,1*2), where spinors Ui and U2 are solutions of (1) 
with the eigenvalues Ai and A2 respectively, 

/l0«l,2 = Ai,2Ml,2 , 

then is a solution of (11). We will choose the spinors Ui and U2 real, which 
is always possible if we suppose Ai and A2 to be real. 

Once the matrix hi is known, the transformation operator L and the trans- 
formed potential Vi are defined up to the matrix A 

L = A[d,-U,U-'], 

Vi = A{Vo + 1^7 - 7^ - -fA-^A^)A-^ . 

The matrix A remains arbitrary provided it satisfies the condition (5). Nev- 
ertheless, remark that the simple gauge transformation (p — A(p reduces it to 
the identity. Therefore, without loss of generality wc can take A = I, where / 
is 2 X 2 unit matrix. Now the operator L has its simplest possible form 

L^d,-U,U-\ (13) 



which is the relativistic extension of the well known Darboux transformation 
operator for the Schrodingcr equation (sec e.g. [25]). For the transformed po- 
tential we get 

Vi^Vq + D, D^VL^-^VL, VL^-U:cU-^. (14) 



We observe here that the transformation operator and the transformed poten- 
tial are completely defined by the matrix-function lA. Therefore, we will call 
it transformation function. 

Using Eq. (11) one can express Ux in terms of U and Vq. This gives us another 
representation for the transformed potential 



V^^-^Voi + D, (15) 
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D = uAu~^ + 7mAm"S = (Ai - A2) A 



d2 —di 



(16) 



where 

di = 1*111*22 + 1*121*21 , d2 = 1*211*22 + 1*ll1*12 , (17) 



A = (detU) ^ and Uij, i,j = 1,2, are the entries of the matrix U. Similar 
transformations apphed to the formula (13) yield 

ip^Lijj^ -f(UAU-^ - E)'ijj . (18) 



This formula is equivalent to (13) provided L acts on an eigenfunction of ho 
with the eigenvalue E. 

Before going further, it is necessary to point out a couple of simple but relevant 
comments. 

Remcirk 1 Given a hamiltonian ho — ^d^ + Vq, if the hermitian potential 
matrix Vo{x) is symmetric, then it can be reduced by a smooth orthogonal 
transformation [22] to the form 

Vo{x) = poix)a3 + qo{x)ai, (19) 



where (Ti,(T3 are Pauli matrices, known as a canonical representation o/ Vq. 
Therefore, in the following we shall use only this canonical representation for 
the Dirac potentials. 

Remcirk 2 Any potential of the form (19) has the properties 

iVo^ = Vo , ^Vo + T/o7 = . (20) 



Note that the operator L has a nontrivial kernel, LU = 0, i.e. kerL is a two 
dimensional linear space spanned by the spinors Ui and U2- This means that 
the set of the spinors Lip, when ip runs through the whole space of solutions of 
the initial equation, may in general not coincide with the whole space of the 
solutions of the transformed equation. In Section 3 we shall show that despite 
this fact a one-to-one correspondence between the spaces of solutions of the 
equations (1) and (2) may be established. 
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3 One-to-one correspondence between the spaces of solutions 

To define the Dirac Hamiltonian and the transformation operator we have not 
used the notion of the Hilbert space and an inner product is not defined in 

the space of spinors. Nevertheless, we need operators adjoint to given ones. 
Therefore in the space of operators we are working we shall define the op- 
eration of conjugation in a formal way. Namely, we will require d'^ = —d^, 
(AB)'^ — i+ = —i (where i is the imaginary unity), and if A is a 

matrix, is its hermitian conjugate in the usual sense (i.e., it is obtained 
from A by transposing and taking the complex conjugation). Moreover, for 
a given invertible matrix U we will assume that the operations of taking the 
inverse and adjoint commute, i.e., 

{u+)-'^{u-y. (21) 

It is easy to see that the Dirac operators ho and hi are self-adjoint. Therefore, 
the adjoint intertwining relation reads 

hoL+ = L+hi . (22) 

This means that the operator L"*" is also a transformation operator and it 
realizes the transformation in the opposite direction, from the solutions of (2) 
to solutions of (1). But as we will show below, it is not an inverse of L. 

It is well known that for any fixed value of the energy the Dirac system (1) has 
two linearly independent solutions. As in the case of the Schrodinger equation, 
if we know one of them, ip with eigenvalue E, the other solution ip with the 
same eigenvalue can be found by a quadrature. In this way we get the spinors 
Ui and U2 with eigenvalues Ai and A2. They do not belong to the space kerL 
and hence vi^2 — Lui^2 are nontrivial solutions of the transformed equation 
with the eigenvalues Ai,2. We will find them next. 

For this purpose we need a counterpart of the Wronskian which, like in the 
case of the Schrodinger equation, may be defined as a function of two solutions 
■0 and ijj corresponding to the same eigenvalue E, which is independent of the 
variable x. Using the adjoint form of the Eq. (1), we easily find 

{^'''y'lp)x = i'li'2x + i'lxi^^ - - i'2xi^l = , (23) 

which means that 

W{'il>, tp) :— 'ip'^^ip — ■0i'02 — ■02'0i — constant . (24) 
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Hence, the function W^'ipjip) we just defined can play the role of the Wron- 
skian. The solutions ip and ip can always be chosen such that the constant 
that appears in (24) is equal to one. 

The Dirac system with a canonical potential Vq = PoCs+Qo^i, written in terms 
of the components ipi and ■02 of the spinor ip — (■0i, ■02)*, has the form: 

'4'2x + PoA + ?o^2 = Eipi , -^13. + qoijji - Pq^2 = Eip2 . (25) 

The same system for the components of the spinor ip = [ipi,ip2y is 

i'2x + PoV'i + ?oV'2 = E^x , -V'lx + qo'tpi - P0^2 = ■ (26) 

If now we eliminate the function from the last equations of (25) and (26) 
we obtain 

ipii^ix - i^ii^ix ^Po + E, (27) 

where we have used formula (24). Prom here, assuming that ipi is not identi- 
cally zero, we easily find 

X 

^i^^PiJ^^^dx. (28) 

xo 



Formula (24) gives us the second component of the spinor ip: 

= — • (29) 

-01 ^1 



If ■ipi{x) = then, for a nontrivial spinor ip = (^i,^/'2)*, 4'2{x) is not identical 
null, and the following alternative formulae should be used: 



xo 



ip2 "02 ip: 



'2 



Now, we can apply (18) to the eigenspinors Ui^2 with the eigenvalues Ai^2 which 
gives us the action of L on the matrix U = {ui, U2). After some simple algebra 
we get 

LU ^ {X2 - Xi)V , (31) 
V=(W+)-^ (32) 
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Hence, the matrix V satisfies the equation hiV = VA, A = diag(Ai, A2). It is 
composed of the spinors Vi and ?,'2. V = {vi,V2), satisfying the Dirac equation 
with the potential Vi, /iifi,2 = Ai^2'*^i,2- Two other cigcnspinors of hi, Vi and 
V2, with respective eigenvalues Ai and A2, /iit'i,2 = Ai^2^'i,2, may be found with 
the help of the same formulas (28)-(30), applied this time to the transformed 
equation. 

Hence, we have established a one-to-one correspondence between the spaces of 
solutions of the equations (1) and (2). For any E ^ Ai, A2 the operators L and 
realize this correspondence; if = Ai^2 the correspondence can be assured 
by the mapping Ui^2 ^ ^1,2, ^1,2 ^ ^1,2, considered as a linear mapping. 

As a final remark of this Section, we notice that L+V = 0, meaning that U E 
kerL+L. In addition, the operator L+ expressed in terms of V reads as L+ = 
—dx — {U^)~^IA^ = —dx + VxV~^ . This means that if we interchange the role 
of the initial and final equations, the function V becomes the transformation 
function for the transformation operator of the type (13), taken with the 
opposite sign and realizing the transformation in the opposite direction, a 
fact that we already mentioned. Moreover, for the transformed equation the 
function V plays the same role that U plays for the initial equation. This 
implies that, according to (31), L^V = (Ai — A2)W. Taking into consideration 
the fact that U G kerL, we find that in one hand V G kerLL"'". On the other 
hand, {U,U} G ker(/io - Ai)(/io - A2), and {V, V} G ker(/ii - Ai)(/ii - A2). 
Hence, we see that 

ker(/io — Ai)(/io — A2) = kerL"'"L = span(W,W) (33) 

and 

ker(/ii - Ai)(/ii - A2) = ker LL+ - span(V, V) . (34) 

From these properties, we can suspect that the operator L^L coincides with 
{ho — Ai)(/io — A2) and the operator LL+ coincides with {hi — Xi){hi — A2). 
These facts will be proved in the next Section. 



4 Factorization property of Dcirboux transformations 

Before establishing the main result of this Section, we need first to prove some 
auxiliary relations. 

Proposition 1 For any real 2x2 matrix A the following formula holds: ■yA + 



9 



Proof. Since the set of matrices {/, 7, cti, 173} is complete in the space of 2 x 2 
matrices, any such a matrix may be presented in the form 

^4 = a/ + 6i(Ti + 627 + ^2,0^ 

and therefore iiA = 2a. The statement of the Proposition follows from the 
fact that the Pauli matrices anticommute. □ 

Lemma 1 The matrix U defined as U := {ui,U2), where Ui and U2 are real 

eigenspinors of Hq (with eigenvalues Ai and A2 respectively) chosen such that 
U is non- degenerate, satisfies the following relationship: 

U := UM~^ - (U^U-y = -(Ai + A2)7 ■ (35) 



Proof. The equation for W (11) together with (21) implies 

{U+)-^U+ = -Vo-f + {U+)-'AU+^ . (36) 

Prom here it follows the relation for U as given in (35) 

U = 'yVo + - -fUAU-^ - {UAU-Y'y > (37) 



which, after taking into account Remark 2, can be rewritten in the equivalent 
form 

u = -Yium-^ - {um-^y^] . (38) 



The statement of the Lemma 1 follows from this equation. Proposition 1, and 
the following equality: tr {UMA~^) = tr A = Ai + A2. □ 

Proposition 2 The following relationships hold for the matrices lA, chosen 
as in Lemma 1, and K, given in (12); 

U^JU'^ = + 7^00. - UA'^U-' , (39) 
A' - (Ai + A2)A = -A1A2/ . (40) 

Proof. The first one is a differential imphcation of the Eq. (11). Indeed, if we 
take the derivative of this equation and replace the first derivative of U using 
the same Eq. (11), we get 

lU^x + VoM + Vo^VoU - Vo^UA = ^VoUA - ^UA^ . (41) 
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Now we use the formulae 

Vo^VoU = -^V^U , VoiUX = -WoUA , (42) 

which are direct imphcations of (20), and rewrite (41) as follows: 

lK.x + VoJA - ^V^U = -tWA^. (43) 



Formula (39) is now evident. 

Formula (40) may be proved straightforwardly, by computing explicitly its left 
hand side using (12). □ 

We establish and prove now the main result of this Section. 

Theorem 1 Darhoux transformation operator L given in (13) and its formally 
adjoint , constructed with the help of the matrixU := {ui,U2), where Ui and 
U2 are real eigenspinors of Hq (with eigenvalues Ai and A2 respectively) chosen 
such thatU is non- degenerate, factorize the following polynomial of the Dirac 
Hamiltonians ho and hi : 



L+L^{ho-XiI){ho-X2l), (44) 
LL+ = (/ii- All) (/11-A2/). (45) 

Proof. Let us consider the superposition L+L. Using (13) and after some simple 
algebraic rearrangements, one gets 

L+L = -dl + Ud^ + U^^ U-^ - UU^ U'K (46) 

We now replace here U by the right hand side of the formula (35) and use 
Eq. (39) obtaining 

L+L ^-dl + Vi + -fVo^ - UA^U-^ 

(47) 

- (Ai + A2)(7fi'. + Vo) -U[A^ - (Ai + A2)A]W-i . 

Formula (44) follows from here, Eq. (40), and the equalities jd^ + Vo — hg and 

-dl + + 7^. = hi 

To prove formula (45), we act with L on the Eq. (44), and take into account 
the intertwining relation (3). As a result, we obtain 

LL+Li) = (/ii - Ai7)(/ii - \2l)L'il) . (48) 
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The last equation means that (45) is vahd for any 99 = Lip and it remains 
to prove it for the spinors which can not be represented in this form. But we 
know that such spinors are eigenspinors of /ii, either with the eigenvalue Ai 
or with A2. For these spinors both right and left hand sides of (45) are zero, 
since the spaces ker LL"*" and ker [hi — Ai/)(/ii — A2/) coincide, as we already 
commented in (34). □ 

As a consequence of this Theorem, taking into account (16) and Remark 2, 
summarized in Eq. (20), the Darboux transformations keep unchanged the 
canonical representation of a potential. 

To conclude this Section, we notice that for the case Ai = — A2 = Aq Eq. (44) 
takes the form 

L+L={h,-\,){ho + \) = {hl-\l), 
a. result that was previously reported in [16] for a scalar potential. 



5 Darboux transformation operators in a Hilbert space 



Let us consider a Dirac Hamiltonian defined in a dense domain D of the 
Hilbert space H = L^(M) ® C^, D <Z H with the inner product 

00 

(^1$):= J il)\x)ip{x)dx . (49) 



Here ip^lx) = {ipKx) , iIj2{x)) . We shall denote the elements of the Hilbert 
space by capital symbols and we shall keep small symbols for the coordinate 
representation of the corresponding spinors. We shall suppose that the differ- 
ential expression + Vq, initially defined in a dense subset D C L'^{R) 
of sufficiently smooth functions, has a self-adjoint extension with domain of 
definition D, which we shall denote by the same symbol Hq. In this case the 
system of eigenspinors G D, /^oI^b) = EI^^^e), is complete in H. This 

means that there exists a measure p{E) such that the following completeness 
condition takes place 

J dp{E)\^E){'^E\^l. (50) 



This equation should be understood in a week sense, i.e., it is equivalent to 

J dp{E){^ri\'^E){'^E\^m) ^Sn,m, n,m = 0,l,..., (51) 
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where |$„), n = 0, 1, ... is an orthonormal basis in H. In coordinate repre- 
sentation the kets \^ e) are reduced to the smooth spinors iJje- Therefore the 
action of L on them is well-defined. Let us consider the following family of 
spinors 

\^e) = Ne^L\^e) , Vl^s) e D, (52) 

where N"^ = (i?— Ai)(i?— A2). It is easy to see that these spinors are elements of 
H if iV|; > 0. Indeed, they are eigenspinors of hi. hi\^E) = E\^e)- Therefore, 
taking into account the factorization property (44), we can define the action 
of L+ on the spinor \^e) 

L+\^e) = AT^^L+LI^b) - Ne\^e) . (53) 



Remark that the right hand side of (49) is nothing but the sum of two integrals 
with respect to Lebesque measure on the real line. Therefore, when L acts on 
eigenspinors of ho and L+ acts on these of hi, the adjoint operation formally 
introduced at the begining of Section 3 coincides with the adjoint with respect 
to the inner product (49). As a result, for all Ai < A2 such that the interval 
(Ai, A2) does not contain spectral points of kg, we have 

($^1$^,) = Ne^{L^e\L^e) = Ne^{^e\L^L^e) = {^eI'^e) . (54) 

This result means that if > A2 or £■ < Ai < A2 is such that E e spec(/io); 
then it also happens that E e spec(/ii). The inverse statement is in general 
not valid. Hence, to find the whole spectrum of hi it remains to analyse the 
values Ai and A2 only. Here, several possibilities can take place. Let us consider 
two real numbers fl and fi, }1 < /i, such that the interval (/i, /x) is a spectral 
gap of ho, for instance the gap between the positive and negative parts of the 
spectrum, if any. Then, the main posibilities are the following: 

(I) If /i is a level of the discrete spectrum, then one can take X2 = (1 < Xi < 
A2, and U2 = i/jf,. 

(A) If, in addition, there exists a linear combination of Vi and Vi, (f — CiVi + 
C2V1 (we recall that Vi (vi) is the first column of the matrix V (V)), such 
that if E H , then the level E = fi will not be present in the spectrum of 
hi, but the new discrete level E = \i will appear. 

(B) If the last condition is not satisfied, the spectrum of hi will coincide with 
the spectrum of Hq, except for one level E — /i, which is missing in the 
spectrum of hi. 

(II) If both Ai and A2 fall inside the interval {fi,ij), /i < Ai < A2 < then the 
following posibilities arise: 

(A) The spectrum of hi totally coincides with the spectrum of ho. 

(B) One discrete level is created at £■ = Ai 
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(C) One discrete level is created ai E = X2. 

(D) Two new levels are created at £^ = Ai and E — X2. 



Since the functions are smooth enough, we know the action of L not only 
on these spinors, but also on any linear combination of them of the form 

I*) = / dp{E)C{E)\^E), (55) 



where C{E) is a finite function over M, i.e., a function with a compact support. 
The set of spinors (55) is a linear space that we will denote in the sequel by 
Cq. We notice that it is dense in H. The image of the space Cq under the 
action of the operator L consists of the spinors 

1$) = J dp{E)NEC{E)\^E). (56) 



This new set will be denoted by Ci. The operator L+, being a linear operator, 
is defined for all |$) e £1. Moreover, the following equality holds 

= (^|L+$), V|^)e/:o, V|<I>)g£i. (57) 



Nevertheless, this does not mean that the operator L+ is adjoint to L with 
respect to the inner product in H. To find such an operator one has to specify 
correctly its domain of definition. We shall not look for this domain. Instead 
we shall give a closed extension L of the operator L and then find its adjoint 
U with respect to the inner product in H. 

For simplicity, let us suppose that hi is completely isospectral with Hq. (If this 
is not the case, the reasonings are similar, but it will be necessary to specify 
the spaces where the operators act on.) In this case the set of functions (52) 
is another basis in the same space H and the operator 

C/ = y dp{E)\^E){^E\ (58) 

realizes an unitary mapping of H onto itself, 

U-'^U^^ J dp{E)\^fE){^E\, (59) 

provided the vectors are orthonormal. Consider now the following oper- 
ators 

L = I dp(E)NE\^E){^E\, (60) 
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I dp{E)NE\^E){^E\, (61) 

It is not difficult to specify their domains of definition. For this purpose we 
introduce the self-adjoint operator go — {ho — Xi){ho — X2), wich has the 
following spectral decomposition 

go^ J dp{E)Nl;\^E){^E\. (62) 

We shall suppose the parameters Ai and A2 to be such that go is positive 
definite. Then, the spectral decomposition of its square root is 

g'J^^ I dp(E)NE\^E){^E\. (63) 



It follows now that 

= ibo^'^ir = / dp{E)Ni,m-fE)\'. (64) 



This means that the domain of definition of L coincides with that of go and 
consists of all ip & H such that the integral in the right-hand side of (64) 
converges. The domain of definition of U coincides with that of the operator 
gl^"^ where gi — (hi — Xi)(hi — A2) is also positive definite. 

From formulae (60) and (61) it follows that the operator L"^ is the adjoint of L 
with respect to the inner product, the domains of definition of L and U being 
well specified. Moreover, Z^^ = L. This imphes [26,27] that the operator L is 
closed. Expressions (60)-(61) give quasispectral representations of the closed 
operators L and U. 

From (60)-(61) it also follows that LipE = LipE = Ne^e and L'I"<^e = L~^(pE = 
Ne^Pe- This means that L is the closed extension of the operator L, and 
is the closed extension of the operator L+, when the domains Cq and Ci are 
taken as their initial domains of definition. 

From the spectral decomposition of the operators g^"^ (63) and 5^^^, 
= j dpNp\(pp){ipp\ , 

one obtains the following representations for L and L+: 
L^Ugl/'^gl/'U, V ^ g^U'^ ^ g^ . 
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Such representations are known as polar decompositions or canonical repre- 
sentations of closed operators (see for example [27,28]). 

Let us consider now bounded operators 

M = y dp{E)NE'\%){^E\ , M^^l dp{E)NE'\^p){^E\ , 

defined in H. It is not difficult to see that both MD and M'^L arc unit 
operators in H. Using the spectral resolutions of the operators Qq and gi 

Qo'^' = I dp{E)N^'\<ilE){^E\ , 9l'^' = I dp{E)N^'\^E){^E\ , 
one derives the polar decompositions of the operators M and M'^: 

It is easily seen that these operators factorize the inverses of g^ and g\. M'^M ~ 
g^\ MM^ = g^\ 

As a final remark of this Section, we would like to notice that the space H 
can be obtained as a closure of the linear space Cq of all finite linear combi- 
nations of the functions = -^I^e) with respect to the norm generated 
by the inner product (49). The set of functions of the form |$) = L\^), 
when 1^) runs through the whole domain of definition of the operator L 
(i.e., the domain of definition of the operator ^/go, D{^fg^)), can not give 
the whole space H. Nevertheless, if one defines a new inner product in £i, 
= (L^'jL^'b) = (^J^ol^fc), l^a.b) e A, \^a,b) e then the clo- 
sure of Ci with respect to the norm generated by this inner product coincides 
with the set |$) = L|*), 1^^) e D{^). This space is embedded in H. 



6 Supersymmetry 

As it has been mentioned in the Introduction, supersymmetric transformations 
in the non-relativistic quantum mechanics are basically equivalent to Darboux 
transformations for the Schrodinger equation. This equivalence is based on 
two properties of the Darboux transformation: (a) the intertwining relations 
and (b) the factorization of the Hamiltonians, similar to that established in 
Section 4. 

In the case of the Dirac equation we started from the intertwining relation 
and proved the factorization properties. Therefore, every property giving rise 
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to the supersymmetry of the Schrodinger equation also takes place for the 
Dirac equation. Moreover, the transformation operators are well-defined in 
the Hilbert space L^(M) (8)C^. Therefore, in order to show the supersymmetric 
character of our approach we can proceed in the same way as in the nonrela- 
tivistic case. 

To begin, let us introduce the following 4x4 matrices 





^ho \ 




(olA 




foo\ 


n = 






Q = 






hi) 




^0 j 







(65) 



It is easily seen that, in one side, the two commutation relations 

[Q,n]^[Q+,n]^o (66) 

are equivalent to the intertwining relations (3) and (22). On the other side, 
the anticommutation relations 

{Q, Q+} = QQ+ + Q+Q = {n- \J){H - X2I) (67) 

are equivalent to the factorizations (44) and (45), while 

Q2^(Q+)2 = 0. (68) 

Remark that relations (66)-(68) are those of a quadratic deformation of the su- 
peralgebra sqm{2), inherent to the usual supersymmetric quantum mechanics 
[8]. This quadratic superalgebra cannot be seen directly from the Dirac equa- 
tion, and therefore we associate it with a hidden supersymmetry. Let us also 
point out that a superalgebra similar to that of (66-68) can also be found in 
the non-relativistic context, when second order Darboux transformations are 
considered [29,30]. 



7 Chains of Dcirboux transformations 

The aim of this Section is to iterate the one-step transformations we have 
considered in previous sections. In order to accomplish this goal, we will rewrite 
first the formulas (13) and (14) in a form wich will be more appropriate for 
this purpose. Remark that the action of the operator L to ^ can be rewritten 
as follows 

= LV' = - W- V = W(W- V)x ■ (69) 
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Now, we will substitute the term in parentheses by the following expression 



W{^,g) 



(70) 



where W{f, g) — p^g — fig2 — f2gi is the Wronskian for the Dirac equation 
that we introduced in Section 4, being / and g the spinors from wich the matrix 
U is composed: W = {f,g), f = (/i,/2)*,5' = (5'i,5'2)*, and hof = Xf,hog = jig, 
i.e., A = diag(A, /i). It can be established by a straightforward calculation that 
the derivative of (70) can be written as 



[W{f.g)Y 



Ws{W{f,g),W{ij,g)) 
\Ws{W{f,g),W{f,i;)) 



(71) 



were Ws{f,g) = fg' — f'g is the usual Wronskian. In the sequel we will use 
both, the prime and the subscript x, for denoting the derivative with respect 
to x; as usual, a multiple derivative of a function / will be denoted as f^'^\ 

After replacing {U~^ip)x for its expression in (71), one gets from (69) the 
desired form for (p: 



W{Lg) 



fi 9i ^1 




/i 9i ^1 


/2 92 




/2 92 P2 


fi g'l ^'i 




n 9'2 i^2 



(72) 



The symbols 



are used to denote determinants. A similar calculation re- 



sults in another representation of the transformed potential matrix Vi given 
in (14): 



14 = 1/0 +[7, A], 



(73) 



18 



where the square brackets represent the commutator. The matrix Di is 



W{f,g) 



fig'i 

f2 92 
f2 92 



h 9i 
f[9[ 

fi 9i 

f^9'2 



\ 



(74) 



In the next subsections we generahzc the formulae (72)-(74) to a chain of n 
consecutive transformations of the same type. 



7.1 Transformation of spinors 



Let us denote now by Li^o the first order operator intertwining ho and hi, 
as given in (13), having the transformation function Ui with the eigenvalue 
Ai = diag(Ai, /xi), Ai ^ /xi. Let U2 be a solution of (1) with A2 = diag(A2, /X2), 
^2 7^ fJ'2, {A2,Ai2} 7^ {Ai,/Ui}. The last condition, together with the non- 
degeneracy of Ui, means the non-degeneracy of the matrix V2 = Li<_o^2- The 
function V2 being an eigenfunction of hi may be chosen as the transformation 
function for the second transformation step. 

The operator carrying out the second transformation will be denoted by ^2^1- 
It intertwines hi and /i2, the new potential given by the formula (14) in which 
the replacements Vq ^ Vi, Vi ^ V2 and — > V2 are done. It follows from 
here that the second order operator 

-^2<-o — L2^iLi^o (75) 



intertwines ho and /i2 and realizes the transformation from ho directly to /12, 
without using an intermediate potential Vi. It is completely defined by two 
transformation functions Ui and U2- It is clear that subsequent iterations give 
us an nth order transformation operator 

Ln<-0 — Ln^n-l " ' " -f'2-(-l-f'l*-0 i (76) 



realizing the transformation between the first ho and the last /i„ elements of 
the chain of Hamiltonians ho, hi, ... , hn. The operator L„<_o is defined by n 
transformation functions Ui, . . . ,Un 

hoUi = UiKi, Ai = diag(Ai, /x,) . (77) 
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We shall derive now a compact formula for the eigenfunctions of the Hamil- 
tonian /i„ (not necessary belonging to the Hilbert space). For the columns of 
the matrix Ui we introduce the notations fi = {fa,fi2y, Qi = (s'ii! fe)*, so 
that Ui = {fi,gi), fi, Qi being two-component vector-columns (spinors), and 
hofi — Xifi, hoQi — HiQi. 

In order to obtain the formulas we are looking for, we need to introduce new 
notations. Let W{fi, gi, . . . , fn, Qn) be a (2n)th order determinant organized as 
follows. The first two rows of this determinant are composed of the components 
of the spinors /i, gi, . . . , /„, gn, the first line from the first components and 
the second line from the second components; any subsequent pairs of rows is 
the derivative of the previous pair. Thus, we can write 



W{h,gi,...Jn,gn) = W{fi,...,gn) = det(aij), i, j = 1,2, . . .2n , (78) 



where 



_ -f('-i) _ ^C-i) 

0'2l-l,2m-l — Jml ) 0'2l,2m-l — Jm2 i 



_ [l-l) _ 
<^2l-l,2m — gml ) C'2l,2m — fl'm2 ; 



(79) 



l,m — 1, . . . ,n. Remember that W{fi, gi, . . . , /„, gn) is the determinant of an 
even order matrix. 



We need also some odd order determinants. If we have the same set of 2n 
spinors fi and gi, i = 1, . . . ,n, plus an additional spinor ip = {ipi, ip2y, then we 
can construct two different determinants of (2n-|- l)th order from the previous 
determinant (78): by adding a column composed of these spinors (using the 
procedure described above) and a row composed of the nth derivative of either 
the upper elements of the spinors /i- (/i- • • . , /„, gn, i^, the lower elements 
of the spinors. In this way, we can get determinants of two kinds 



Wi{fi, gi,..., fn, gn, ip)^det{bi,j) , (80) 
W2{fi,gi,...,fn,gn,tp)^det{cij), i,j = 1,2, ...2n + l, (81) 
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(82) 



where bij — Cij — when — 2n, and 

b2l-l,2n+l — C2l-l,2n+l —'4^1 ^\ I — 1, ... ,71] 

h2,2n^.=C2,2n^.=^t'\ l = l,...,n; 

1 r(n) f{n) -I 

02n+l,2m-l — /ml > C2n+l,2m-l — Jm2 ) — 1, . . . , n; 

1 (n) (n) 

02n+l,2m — 9ml) C2n+l,2m-l — Qna-i — i, ■ ■ ■ , 

^2n+l,2n+l ="01^ C2n+l,2n+l — "^2 ^ ■ 



Now, in order to prove the main result of this Section, we need the following 
auxiliary statements. 

Lemma 2 The following equation takes place: 



Ws{W{h, g,,..., /„, g^), W{f,, ^i, . . . , /„, ^)) = 

= Wi{fi, gi,..., fn-l, gn-l, /n)M^2(/l, ■ ■ ■ , fn, Qn, "0) (§3) 
-W2{fl, gi,..., fn-l, Qn-l, fn)Wi(fi, gi, . . . , fn, gn, "0) ■ 



Proof. To prove the lemma we need first the following identity, which can be 
checked in [19]: 

W^eW' - W'Wl, = W'^'We-,, (84) 



where 

:= Wi,2,...,e-i,r,i,2,...,e-i,k, 3, k > £, 

^h,i2,--;ik,;ji,j2,--;jk ■~ ^^^{yimjn), TTljU = 1, . . . ,k, 

and Ui^j^ are arbitrary. 

Note first that the derivative of a determinant such as the one defined in (78) 
is the sum of two determinants of a similar structure. The difference is only 
either in the last row or in the next-to-last row. In the first case we obtain the 
determinant denoted by Ri{fi, gi, . . . , fn, gn) in which the last row is replaced 
by the nth derivative of the second elements of the spinors /i, gi, . . . , /„, gn 
and in the second case we get the determinant denoted by i?2(/i, S*!, • • • , fn, Qn) 
in which the next-to-last row is replaced by the nth derivative of the first 
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elements of the same spinors. To be more precise, we have the following 

[W{fu Qn)]' = Rlifl, gi,---, fn, 9n) + ^2(/l, ^1, • • • , /n, Qn) , (85) 

where 



Ri(fi,gi,---,fn, Qn) = det(r- J-) , 

i?2(/i, • • • , /n, 9n) = det(r? •) , i,j = 1,2,.. .2n, 



and 



= ^Ij = ^^.i' ^ = 1, • • • , 2n - 2, J = 1, . . . , 2n; 

1 An) I (n) -> 

^2n-l,2m-l ~ /ml ) ^2n-l,2m ~ fl'ml ) — i, . . . , n; 

2 2 C*^) 1 

^2n,2m-l ~ Im2 : ^2n,2m ~ 9m2i TH — 1, . . . ,n; 

'^2n,j — ^2n,3, ^In-lJ — «2n-lJ, J = 1, ■ ■ ■ , 2n. 



(86) 

(87) 



The functions defined in (79). The Lemma follows now from (84) and 

(85). ' □ 

Proposition 3 In what follows, we will need also the following identity 



Wiifi, gi,..., fn-l, 9n-l, fn)W2{fl, Ql, ■ ■ ■ , fn-1, Qn-l, Qn) " 

-Wiifi, 91,..., fn-l, 9n-l, 9n)W2{fl,9l, fn-1, 9n-l, fn) (89) 

= W(fi, 9l,---, fn-l, 9n-l) W{fi, 9l,..., fn, 9n) ■ 

which is a direct consequence of {84). 

Now we formulate and prove the main result of this Section. 

Theorem 2 The action of the operator Ln<-o on a function ip = {'4^i,'4^2f 
reads as follows 



I 



W{fi,gi,...,fn,9n) 



Wi{fi,gi,. . . ,fn,9n,ij) 



^W2{fl,gi,...,fn,gn,'^) J 



(90) 



where Wi and are defined by (80)-(82). 
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Proof. To prove the theorem we use the perfect induction method. Let us 
suppose that the action of the operator L^-i^-o on a function -0 have the form 
(90), with the replacement n — > n — 1, i.e., 



^ W2(fl, gi,..., fn-l, Qn-l, -0) y 



(91) 



Then, according to (76) and (69), we have 



(92) 



where 



Un — Ln-l^oUn , Un — {fn, 9n) 



(93) 



Here, the spinors and (?„ should be calculated by the same formula (91), 
where ip has to be replaced by /„ and gn, respectively. Using Eq. (70) we find 



1 



W{fn,gn) 



^W{iP,gy 



(94) 



After calculating those Wronskians, and using (91) and (89), we get from (94) 
the equation 



W{fi,gi,...Jn:9n) 



-W{fi,gi,...,gn,^) 



^ W{fi,gi,...,fn,ip) J 



(95) 



The derivative of this function reads 

/ 



1 



[W{h,...,gn)Y 



Ws{W{h,...,gn),W{h,...,gn,^l^)) 
[ Ws{W{h,...,gn),W{h,...,U,i^))J 



The statement of the Theorem is a direct implication of last equation, together 
with equations (83), and (92). □ 
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1.2 Transformation of the potential 



In the last part of this Section we need to introduce the following nota- 
tion: let us denote by Qi{fi,gi, • • • , fn, Qn) the determinant constructed from 
9i7 ■ ■ ■ 7 fn, Qn) by the replacement of the (2n)th line with the nth deriva- 
tives of the first elements of the spinors /i, gi, .... /„, and let us denote 
by Q2{fi,gi, • • • , /n, Qn) the determinant obtained from W{fi, gi, . . . , fn, gn) 
by the replacement of (2n — l)th line with the nth derivatives of of the second 
elements of the same spinors. Thus, we have 



Qi (/i, ^1, ■ ■ ■ , fn, gn) = det(g- , 
(32(/i, gi,---, fn, gn) = dct(gjj) , i,j = 1,2,.. .2n, 
where 

^Ij = = ^ = 1, • • • , 2n - 2, j = l,...,2n; 

1 fin) I (n) 

?2n,2m-l ~ J ml , Q.2n,2m ~ 9ml, TTl — 1, . . . ,71] 

2 r{n) 2 (") i 

l2n-l,2m-l ~ Jm2 , Q2n-l,2m ~ 9m2, ^ — i-, ■ ■ ■ ,71; 



(96) 

(97) 



(98) 



and the functions defined in (79). 

Theorem 3 The potential resulting from a chain of Darboux transformations 
(76) has the form 



K = ^0 + [7, Dn] , 



(99) 



where 



W{fi,...,gn) 



Ri{fi, gi, ■ ■ ■, fn, gn) Qiifi, gi,---, fn, gn, 



\Q2{fl,g\, ■ ■ ■ , fn,gn) R2{fi,gi, ■ ■ ■ , fn,gn) I 



,(100) 



being R\ and R2 defined by (86)-(88); and Qi and Q2 defined by (96)-(98). 

Proof. To prove this Theorem we use again the perfect induction method. 
Therefore, let us suppose that the formula (99) is vahd for the potential 
Then, Eq. (14) implies 



K = K-i + [7, {K').K'] = Vo + [7, Dn-i + {U-%U-'] . 



(101) 
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Here the matrix Un defined in (93) has the form 



Wr. 



Wi{fi,gi,..., fn-l, 9n-l, fn) ^1, • 

M^2(/l, gi,---, fn-l, 9n-l, fn) W^2(/l, 9l, ■ 

The derivative of this function reads 



V 



; fn-1, 9n-l) 9' 
) /n-1) 9n-l, 9n) j 



(102) 



.(103) 



[W{fi,gi,...,fn-u9n-l)]' 

W{fi,gi,...,fn-l,9n-l) 
1 ,^ . 

^W{h,g„...,U-i,gn-iy''''' 



and hence 



(104) 



(105) 



/7l \ 7~J-l- [W{fl,gi,...,fn-l,gn-l)]' j 

\L4n)x l^n — T7777 7 \~ 

W{fi,gi,...,fn-l,gn-l) 

W{fi,gi,...Jn-l,gn-l) 



After calculating the derivatives and using the same technique as while proving 
the Lemma 2 we obtain 



I 



1 



Rlifl, gi,---, fn-l, gn-l) Ql{fl,gi, ■ ■ ■ , fn-l, gn-i) 



+ 



W{h,.. 


,gn- 




[W{fi, . . 


■ ,gn 


-.)!' 


w{h, . . 
1 


■ ,gn 


-l) 




,gn- 


4- 


1 




(r 



W{f,,...,gn) 



^Q2{fi, gi,---, fn, gn) R2{fi, gi,---, fn, gn) I 



(106) 
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The sum of the first and the third terms in this expression gives 

gU---, fn-l, gn^l) + i?2(/l, gi,---, fn-1, gn-l)]/W{fi, gn-l) I. 

According to (85), this is equal to the derivative of W{fi, . . . , gn-i), and hence 
these items cancel out when added to the second term of (106). Therefore, we 
get exactly the formula (100). □ 

As a final remark of this Section, we notice that the formulae (90) and (99)- 
(100) can be considered as relativistic analogs of the Crum determinant for- 
mulae [24] , a result which is very well-known in the non- relativistic case. More- 
over, since they present the final result of the action of a chain of first order 
transformations, the operator L„<_o and its formally adjoint L'^^q satisfy the 
following factorization properties: 

n 

L^^oLn^o = n (^0 - A,) {ho - fij) , (107) 

n 

K^oL^^o = IK^n - \j){K - ■ (108) 
This result comes out directly from Theorem 1. 



8 Pseudoscalar potentials 



8. 1 Darboux transformation for a pseudoscalar potential 



A general pseudoscalar potential is defined only by one function qo{x), x G 



' m qo[x] 



^qo{x) -m y 



(109) 



where m is the mass of the particle. The Dirac system for the components of 
the spinor = {ipi,ip2y is 



+ = (-E - m)V'i ■ 



(110) 

(111) 



We would like to notice the following property of this system, that we will 
use in the sequel: if one of the components of the given spinor is zero, then 
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the other is also zero for all values of E, except if = ±m. When E = m 
the system has a solution of the form -0 = (-01, 0)*, and when E — —m the 
solution is = (0, ■02)*- 

In general, after applying the Darboux transformation to a pseudoscalar po- 
tential we get a potential which is not pseudoscalar anymore. Here we shall 
formulate additional conditions for Darboux transformations to keep pseu- 
doscalarity of a potential. 

It is easy to see that if one of the elements of a transformation function is zero, 
then the value di defined by (17) is constant, di = ±(Ai — A2). This means 
that after such a transformation the new potential is also pseudoscalar, the 
role of the mass being played by di — m. 

Let us take one of the component of the spinor ui equal to zero, for instance 
U21 — 0. This is possible for Ai — m. In this case 



U 



^ Uii U12 



U22 



(112) 



detZY = uiiU22-i and the potential Vi, given in (15), takes the form 



Vi = -A2C73 + 



.U12 

(A2 - m) go 

U22 



(71. (113) 



Now, from the Dirac system we can find the function 1^12 

U12 = 1-^ (^22 + ^0^*22) , (114) 

Ao — 771 



and rewrite Eq. (113) as follows 

Vi = -A2(73 + giCTi , (115) 

where 

gi = ^ = (1117.22)'. (116) 
1*22 

In the transformed Dirac system the role of the mass is played by — A2. We 
would like also to mention a relationship existing between the nonzero com- 
ponent of the spinor ui and the potential qo, which immediately follows from 
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the initial Dirac system (110): 



u 



11 



Uii 



(117) 



Let us find now solutions of the transformed equation. We first calculate the 
product UxU~^: 



(Iniiii)' ' 



V 







UuU22 

(In 1*22)' 



Then, we simplify this expression with the help of equations (114), (116) and 
(117): 

(Infill)' -(A2 + m) 



^ /^go-(A2 + m)^ 



V 



[0 q, ) 



(lnM22)' y 

Finally, using (13), we find the action of the operator L on solutions of the 
initial equation: 



( {\2-E)i,2 ^ 



-02 - (ln'U22)''02 



(118) 



We observe that the lower component of this spinor is defined just by the same 
expression that appears in the Darboux transformation for the Schrodinger 
equation with the transformation function ^22 (see e.g. [25]). We also notice 
that the upper component of the spinor (118) differs from ■ip2 only by a constant 
factor. This means that they should satisfy the same equation. Later, we shall 
show that this is really the case. 

For the case uu = and A = —m, similar calculations give us the following 
result: 



Vi = A2(T3 + qiai , 
gi = -(lnMi2)', 



'^Vi-(ln^i2)>i^ 



(p-. 



(119) 
(120) 

(121) 



\^ (A2 + E)ilji J 
Here the mass in the transformed Dirac system is equal to A2. 
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As a conclusion, we have shown that for a pscudoscalar potential both the 
transformation operator and the transformed potential are expressed in terms 
of just one function. 



8.2 Interrelation between Darboux transformations for the Dirac and Schrddinger 
equations 



It is well known (see e.g. [31]) that, when using a pseudoscalar potential, the 
Dirac system may be reduced to the following two independent Schrodinger 
equations 

+ Vi^eV'i, (122) 

-V'^' + C/(-V2 = £V'2, (123) 

where 

C/f ) = ql ± (124) 

and e — E^—m^ . The system (122)-(123) is precisely a pair of super symmetric 
Schrodinger equations (see e.g. [32]), one equation being the SUSY partner of 
the other. Similarly, one of the Schrodinger Hamiltonians 

H^^^ = -d,. + Ut^ (125) 

is the SUSY partner of the other. Moreover, according to (117) one has 

Ut^^ut^-2{]nun)", (126) 

where Uu is an eigenfunction of H^'^^ which is everywhere non-vanishing. The 
transformed Dirac equation is also an equation with a pseudoscalar potential. 
Therefore, it can also be reduced to the pair of supersymmetric Schrodinger 
equations 

-(^;' + C/^+Vi=£i<^i, (127) 
-^'^ + Ut'>^2 = ei^2, (128) 

where 

C/W = qf ± q[ , (129) 
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El = E'^ — \\, and qi is given by (116). We see that to this system corresponds 
an energy different from the one that appears in equations (122)-(123). To 
compare this system with (127)-(128), we have to displace the energy Si. For 
this purpose we add to the left and right hand sides of the equations (127) 
and (128) the terms (A2 — m?)ipi and (A2 — rn?)ip2^ respectively. This leads 



to shifting the potentials f/} 



(±) 



(±) 



(±) 



m 



Now, taking into 



account that U12 and U22 are solutions of the system (122)-(123) with the 
potentials (124) and the expression (116) for the potential qi, we get 



Ut^ = -2{\nu22)" + ql ~q'o = f/^ - 2{\nu22)" 



(+) 



(-) 



(130) 
(131) 



We observe that the potentials Ui^'' and f/g coincide and, as it has been 
mentioned in the preceding Section, the functions ipi and 1IJ2 satisfy the same 
equation. 

Using the expression (117) for the potential go we obtain the potential differ- 
ences 



r(+) 



-2(lnM22)", 
-2(lnMii)", 



(132) 
(133) 



that agree with (126). Hence, we can obtain the potential u[ ^ by two different, 
but equivalent, ways: 

• The first possibility is to start with the initial Dirac system, realize the 
Darboux transformation with the transformation function U given by (112), 
and then split the resulting Dirac equation into a system of two Schrodinger 
equations, related to each other by another Darboux transformation. This 
corresponds to the path qo ^ qi in the diagram of Figure 1. 



Mil 




Fig. 1. Diagram showing the connections be- 
tween the Dirac system, its transformed, and the 
associated SUSY Schrodinger equations. 
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• The second possibility is to split the Dirac equation into two Schrodinger 
equations with potentials U^, and then realize a chain of two transforma- 
tions at the level of the Schrodinger equations, starting with the potential 
Uq] for the first transformation we use uu as the transformation function, 
and get the potential Uq; finally transforming this potential with the help 
of the transformation function U22, wc obtain the same potential ui \ This 
path corresponds to qo Uq^^ Uq~^ Ui~'^ in the diagram. 

Other possible path is Qq — > Uq^^ — > Ui^^ — > Ui~\ which is completely equiva- 
lent to the previous ones, as it is clear from the preceding discussion. Moreover, 
as it follows from the scheme of Figure 1, the following proposition holds. 

Proposition 4 A Darboux transformation for the Dirac equation with a pseu- 
doscalar potential (109), with transformation function (112), is equivalent to a 
two-step SUSY transformation for the Schrodinger equation with the potential 
Uq given by (124) with transformation functions Uu andu22- 

For the case Ai = — m, Uu = 0, similar calculations give us 



We conclude this Section by stressing that Darboux transformation for the 
Dirac equation induces Darboux transformations for the associated Schrodinger 
equations. 



9 Scalar potentials 

9. 1 Darboux transformation for a scalar potential 

In this Section, let us consider that the initial potential Vq has the scalar form 
Vq{x) ^ PQ{x)ari , po^m + So{x), (136) 

where m is the mass of a particle and the function Sq{x) is supposed to be 
known. Note first that if the function ip is a solution of the Dirac equation 
with the potential (136) and the energy E, then the function aiip is also a 
solution, but with the energy —E. In particular, this means that the spectrum 
of the Dirac equation is symmetric with respect to = 0. 



c7|+) = -2(lniii2)" + ?o' + ?o- 



(134) 



(135) 
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In the literature another representation of scalar potentials is frequently used: 

Voix) = {m + Soix))ai . (137) 



Both representations (136) and (137) are related by a unitary transformation 
U: 

Vo = U-'VoU, iP = Uilj, (138) 



where 

U^{l + ^)/V2. (139) 

Observe that the scalar potential (137) can be considered as a pseudoscalar 
potential (109) when the value of the mass is equal to zero. Nevertheless, it 
has some new properties to be considered below. 

The Darboux transformation, in its general form, does not preserve the scalar 
character of a potential. Therefore, it is necessary to select those transforma- 
tions such that they will preserve the scalar form of a potential. As it follows 
from (15)-(16), the transformed potential remains to be scalar if d2{x) = 0. 
From (17) one easily sees that to satisfy this condition it is sufficiently to 
construct the transformation function U — (mi,M2) from the spinors ui — 
{uii,U2iy and U2 — criUi with eigenvalues Ai = A^°^ and A2 = — A^°\ respec- 
tively. 

In order to look for new properties of the Darboux transformation for the 
Dirac equation with a scalar potential, it is more convenient to consider the 

scalar potentials written in the form (137). Under the transformation (139) 
the potential Vq goes into Vq — {111 + SQ)ai, and the transformation function 
U into U, given by 



\U21 U21 J 



(140) 



where un = un — U2i, U21 = Uu + U2i- For calculating the new potential, 
we apply (15)-(17), and for getting the transformation operator we use (13), 
where we have to make the replacement — > ZV. The matrix UxU~^ is now 
diagonal: 



^(InSii)' ^ 
^ {\nu21)' J 



(141) 
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Prom Eq. (14) we find that the transformed potential is 



(142) 



where 



Si = So + (InSai)' - (Inwn)'. 



(143) 



Solutions of the Dirac equations with the potential (142)-(143) are found by 
the action of the operator L = d^—Ux on solutions of the initial equation: 



i)[ - (InSii)'^i 

V'2 - (ln{t2l)''02 



(144) 



Remark that if in the positive part of the spectrum of the initial Hamiltonian 
there exists a ground state level E — Eq with wave function iI)q{x), then in 
the negative part of it there exists a highest energy level E = —Eq with wave 
function aiipQ^x). The choice of the spinors Ui = ipQ and U2 = (Ji'ipo generates 
a potential Vi with the same spectrum as Vq, except for the levels ±Eq. 

The use of formula (140) gives us the matrix solution of the transformed Dirac 
equation with the matrix eigenvalue A 



V= (W^ 



u 



11) 



[U 



\ 



11/ 



{V1,V2) 



(145) 



We observe here that the spinors vi and V2 are either square integrable or 
non-integrable simultaneously. In the first case the levels E — ±A*^°^ appear 
in the spectrum of hi. Hence, the Darboux transformation may create the 
energy levels only by pairs symmetrically disposed with respect to = 0. 
This agrees with the fact that it produces a scalar potential which may have 
only a symmetrical spectrum. 



9.2 Interrelation with the Schrddinger equation 



It is well-known (see e.g. [33]) that the Dirac system with the scalar potential 
(137) may be reduced to the following supersymmetric pair of the Schrodinger 
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equations: 

-iP'l + UM^i = si^i , 

(146) 

Since the transformed Dirac system also corresponds to a scalar potential, it 
may be associated with a couple of similar equations 

-(p'{ + = £1^1 , 

(147) 

where 

[/W = (m + 5i)2±5i. (148) 

Taking into account the equation for S\ (143), after some algebra we get from 
(148) the potentials 

[/(+) = (m + S^f + 5^ - 2(lnSn)" , (149) 
= (m + ^o)' -S'^- 2(lnM2i)" • (150) 

Hence, we also conclude that the potentials u'^^ are SUSY partners of the po- 
tentials f/g^'' , and Darboux transformation for the Dirac equation induces Dar- 
boux transformations of corresponding supersymmetric pair of initial Schrodinger 
equations. 



10 Illustrative examples 



In this Section we will show how the technique we have developed for Dirac 
systems is applied to some interesting examples of pseudoscalar and scalar 
potentials, as well as to spherically symmetric potentials. 



10.1 Pseudoscalar potentials 



First, we will analyse some examples of Darboux transformation applied to 
transparent potentials, and also to the relativistic harmonic oscillator. 
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10.1.1 Transparent potentials 



In this paper we do not analyse the changes in transmission and reflection 
coefficients produced by Darboux transformations. Nevertheless, we can eas- 
ily notice that if the initial potential is transparent, i.e., if it produces the 
zero reflection coefficient for an incident particle, or if the potential does not 
change the asymptotic form of a continuous spectrum cigcnfunction, then 
the transformed potential keeps this property unchanged. In particular, this 
means that starting with the free particle Dirac equation we shall get only 
transparent potentials. Hence, in this Section let us consider the potential 



(151) 



which is a particular case of (109) with go = 0. 

Example 1 Let us take as transformation function the following matrix 



U 



sinh kx 



cosh kx 



k — \lw?- — 



(152) 



It corresponds to Ai = m, A2 = £ < m. The transformed potential is obtained 
by Eqs. (115) and (116), that generate the well-known one-soliton potential 
(see e.g. [34,35]): 



Vi = —£(73 -|- k tanh kx ai 



(153) 



From Eq. (152) we find the matrix solution of the transformed Dirac equation 
for the particular value E — A — diag(m, e): 



I 



V = (w+)-^ 







tanh /ex sechA;x 



(154) 



We conclude from (154) that the potential (153) has one discrete level at _E = 
e. The solutions of the Dirac equation with the potential (153) at E 7^ £:,m 
may be found by applying the operator L to solutions -0 of the free particle 
equation: 

Lp = Li) = i)^-U^U-^-^. (155) 



The potential (153) may be considered as the initial potential for the next 
transformation step. To carry it out, we need a spinor solution of the Dirac 
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equation with either upper or lower component equal to zero. If we take as a 
solution of the free particle Dirac equation 



^ — r— sinh kx ^ 
cosh kx 



(156) 



then, from (155) we find the following spinor solution of the Dirac equation 
for the potential (153): 



^ -2£coshA;x^ 



(157) 



corresponding to the eigenvalue Ai = —e. Another solution of the Dirac equa- 
tion for this potential (153) may be found with the function 



1p2 



( 



\ 



e\—m 



-k\x 



I 



(158) 



in (155), which gives us 

' (£-£i)e-'=i^ ^ 
^ (/ci + A; tanh kx)e~^^^ j 



^2 



ki 



m? — 



(159) 



This spinor has the eigenvalue A2 = ^i- From the spinors Ui = ipi and U2 = v^2 
we construct the matrix solution for the potential (153), U = (mi,M2), with 
the matrix eigenvalue A = diag(— £, £1). 

Example 2 When the above matrix solution U is taken as the transformation 
function for the second transformation step, it produces a two-soliton potential 
of the form: 



V2 = -£1(73 



ki + k tanh kx 



k tanh kx] ui . 



(160) 



The choice /ci > /c > assures the regular behaviour of this potential V2, 
which keeps the discrete level E = e unchanged and has an additional level at 
E — —e. The last statement can be easily seen from the matrix solution for 
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the potential (160) with matrix eigenvalue A = diag(— £, £i): 

( ^— ^ 

2f nosh fc.T 



\ 2£coshfea;(A;i+fetanhfea;) fei+fetanhfea; / 



The first column of this matrix is a square integrable spinor with eigenvalue 
E^-e. 



Example 3 In (155) let us choose now 



1p2 



^ — ^ cosh/cix^ 

£i— m 



V 



sinh kix 



rn? — e\. 



(161) 



This gives us 
I 



<P2 



{e — Si) srnhkix 
ki cosh kix — k tanh kx sinh kix 



, k = \Jm^ — , 



(162) 



The use of this spinor as the second component of the transformation function, 
U2 = (p2, together with the previously found Ui as the first component, 14 — 
{ui,U2), produces the following three-soliton potential: 



Vo 



kl-k'' 



ki coth kix — k tanh kx 



k tanh kx] ai , 



(163) 



which is regular provided ki > k > 0, and has three discrete levels: E = ±s 
and E — Ei. 

It is important to stress that similar potentials have been found recently by 
other means [34] . In contradistinction to these authors, we are able to indicate 
precisely the position of the discrete levels. Moreover, corresponding wave 
functions are easily obtained from columns of the matrix- function (U~^)~^ 
In the next example we give more general transparent potential with three 
discrete levels. 



Example 4 We can take the following spinor as a solution of the free particle 
equation in order to get a solution of the Dirac equation with the one-soliton 
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potential, which is: 



^ ei-m (si^h kix + B cosh kix) ^ 



y (cosh kix + B sinh kix) j 



(164) 



where B is an arbitrary constant. From (155), we obtain 

^ (£-£l)Q(x) ^ 

^2 = 

^ fci(sinh kix + B cosh fcix) — kQ{x) tanh fca;) y 



(165) 



where Q{x) = cosh /cix + B sinh fcix. When ip2 is used as the spinor U2 for the 
second transformation step, it produces the potential 



Vo 



+ 



{kl - k^){l + BtcLnhkix) 



A;i(tanh kix + B) — k tanh kx{l + B tanh kix) 



(166) 
k tanh kx \ ai. 



It is not difficult to prove that it is regular provided ki> k > Q and -B > 1 , and 
also that it has three discrete levels: E = ±e and E = si. A typical behavior 
of this potential is shown in Figure 2, where the term between parentheses 
in Eq. (166) is plotted. One can notice that closer is the parameter B to the 
value 1, wider is the potential barrier. A potential of this type is not known 




Fig. 2. Three-level transparent pseudoscalar po- 
tentials with m = 1 and e = 0.5. Curve 1 is 
obtained for B = 1.5 and e = 0.45 Curve 2 is 
obtained for B = 1.000005 and e = 0.3. 

in the available literature. 

To finish this subsection, we would like to remark that our method can produce 
transparent potentials of a more general form, which are a superposition of 
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scalar and pseudoscalar potentials. We illustrate this fact in he next example 

Example 5 As a final case of transparent potentials, let us consider as trans- 
formation function the following matrix solution of the free particle Dirac 
equation: 



' r- sinh kx ^ 



U 



cosh kx 



-,kx 



, Ai = — £, X2 — e, k = \/m^ — z^. (167) 



It generates a completely new potential of the form 



m — 



ee 



2e 

2kx 



+ m 



(J3 + 2ke 



sinh kx — cosh kx 



^gfea; _|_ 



-kx 



(7i . 



(168) 



From the analysis of the function 



(W+) 



k 



-kcoshkx \ 



g^kx _|_ lYiQ-kx se'^kx _|_ 

171 + e {e — m) sinh kx 
we conclude that Vi has one discrete level E = —s. 



10.1.2 Darhoux transformation for the Dirac oscillator 

As it is well known, in the literature there arc several possible candidates to 
be considered as the relativistic Dirac oscillator. We choose the one which is 
described by the Hamiltonian [36] 



X 



(169) 



Its discrete spectrum consists of a positive series 



(170) 



and a negative series 

E^-^ ^ -{rri^ + ny/\ n = 0, 1, . . . 



(171) 
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Wc will use this model to illustrate the aplications of the Darboux transfor- 
mation for the Dirac equation in three examples. 



Example 6 We will use the spinors 



V / 



(172) 



and 



U2 = 



He.(ix)' 
^ e"^"/"^ Hen_i(ix) ^ 



(173) 



for constructing the transformation function U = (ui,M2)- The spinor (172) 
is a solution of the Dirac equation with the potential (169) for Ai = m; the 
same is true for the spinor (173) with eigenvalue A2 = £n = (^^^^ ~ n)^/^. Here 



(174) 



where Hn{z) are Hermite polynomials. Using the formulas (115) and (116) we 
obtain the transformed potential 



X 



1) 



(175) 



where Kn{x) — (— i)" He„(ix). Observe that the functions Kn{x) are real and 
nodelcss for even values of n; for odd values of n they have only one node at 
a: = [29]. Therefore, the potentials (175) are real and regular when n takes 
odd values. The analysis of the function (W+)~^ shows that this potential has 
two additional discrete levels E — m and E = Sn with respect to the initial 
harmonic oscillator potential (169). 

Example 7 Let us keep the spinor U2 as in the previous example and let 
us take ui = (0, /^)*, which corresponds to Ai = —m. Using the same 
procedure of the previous example, we get the potential 



X Kn-l{x)\ 



(176) 



In contrast to the previous example, Ai belongs now to the discrete spectrum 
of the initial potential. Therefore, the level E — —m is deleted from the 



40 



spectrum of the Hamiltonian (169) and the new level = e„ is added. The 
potential (176) is everywhere regular for even values of n. 

Example 8 Finally, let us take now the spinor ui as in Example 1 and choose 
the spinor U2 as follows: 



U2 



( 



(177) 



which corresponds to A2 = [rn? — 1)-^/^. Here, Q{x) = Yvr/2 {B + erf(a;/v^)), 
with \B\ > 1. We get for the transformed potential 

Vi = -Asaa + (| + e'^'/^ Q-\x)^ a, . (178) 



A simple analysis shows that this potential has two additional discrete levels 
E = m and E = [m^ — 1)^/^ with respect to the initial potential (169). Two 
typical representatives of this potential are displayed in Figure 3, where the 
term between parentheses in Eq. (178) is plotted. It is clearly seen that closer 
is the parameter B to the value 1, larger is the perturbation of the initial 
potential. 



1 


2 


q ^ / 




X 


-4 


/ ~v ^ 

-1 

■^'^ -2 


2 4 



Fig. 3. Harmonic oscillator potential 
q^ix) = x/2 (curve 3) with m = 1, and 
two of its Darboux transformed partners qi{x). 
Curve 1 is obtained for B = 1.0002, and curve 2 
for B = 1.2. 



10.2 Scalar potentials 

In this subsection we will analise some scalar potentials of the form (137) 

Vq{x) = po{x)ai = (m + So{x))ai . 
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In the first instance, we will consider the simplest case, corresponding to 
So{x) = 0. 



10.2.1 Transparent potentials 
Example 9 The spinor 



Ui 



m—k kx I m+k „—kx 
e -r ^ e 



k = Vm? — A^, 



(179) 



is a solution of the Dirac equation with potential Vq = mai for E = \ < m. 
Choosing the transformation function in the form lA = {ui,U2), U2 = —a^Ui, 
and using Eq. (143), we get the transformed potential 



2k' 



m + X cosh(2/i;x + 2a) ' 



:i8o) 



where e^" — yj(m — k)/{m + k). This is a transparent potential, which has 
been previously found in [33]. It is easy to see that the two spinors coming 
from the matrix {U^)~^ are square integrable. This means that the potential 
(180) has two discrete levels: E — ±A. 



Example 10 Consider now the spinor 



/ 



sinh kix 



^^sinh(A;ia; + 2ai) j 



Im — ki 
m + /ci ' 



(181) 



satisfying the free Dirac equation for £■ = Ai < m. After acting on it with the 
transformation operator of the Example 1, we get the spinor 



Wi 



ki cosh kix — k tanh kx sinh kix 
^ki cosh(A;ia; + 2ai) — ktanhkx sinh(A;ia; + 2ai) ^ 



(182) 



which is a solution of the Dirac equation with the potential (180). If we choose 
the spinors Wi and = —CaWi for the next transformation step, we obtain a 
two-step potential 
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Q _ ^1 ^ 

ki coth(A;ix + 2ai) — k tanh(A;x + 2a) 

. 183 

ki coth kix — k tanh kx 

For < k < ki < m this is a regular transparent potential with four discrete 
levels E = ±A and E = ±Ai. Figure 4 shows the typical shape of such 
potentials. From this figure it is clearly seen that closer are the discrete levels 
of the potential, more distant from each other are the two potential wells. 




Fig. 4. Four-level scalar potentials with m = 1 
and A = 0.6. Curve 1 is obtained for Ai = 0.58, 
and curve 2 for Ai = 0.2. 



10.2.2 Scalar Coulomb potential 

We will analise now radial potentials of the form 

Vo = [m ^ cTi, a > 0, 



(184) 



which find applications in modelling inter-quark interactions [37]. In this case, 
the discrete spectrum of the Dirac Hamiltonian consists of a positive and a 
negative series [38] 



Er, 



\ 



1 - 



a" 



'n + aY ' 



n= 1,2, 



(185) 



plus the zero energy level E = Q [39]. The eigenf unctions of the discrete 
spectrum (not normalized here) have the form 



(2a)„ ^ -^n \^-^) 



(186) 
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where = y fn'^ — and x = Enr, n = 1,2, . . . . 

Example 11 In order to apply the two-step Darboux transformation de- 
scribed in Section 7, let us take the spinors U = {ui,U2), with ui = Tp^ and 
U2 = -o-s-^fc, and V = (^1,-^2), with vi = ^^+1 and V2 = -a^-i/jk+i- For the 
transformed potential we obtain 

Vi'\x) = {m + Si'\x))a^, = -- + (In g^W)' - (In gi(r))', (187) 



where 



Qi{r) 
Q2{r) 



kLl-^\2e,^,r)Ll--\2e,T) {k + l)Ll-_\\2e,T)Ll%,\2e,^^ 



r] 



2a + k + l 2a + k ' 

i2e,^ir)L'r' {2e,r) Ll^_f {2e,r) L^^^ {2e,+,r) 
{a + k + iy {a + kf 



The spectrum of V2^\x) differs from the spectrum of the initial Coulomb 
potential by the absence of the levels E = Ey. and E = E^^i. We would like to 
remark that after the first transformation, either with the spinor Ui or with €,2, 
we get potentials with singularities, but the second transformation removes 
all singularities and the potential (187) is regular for r > 0. The simplest 
particular case corresponds to A; = 1: 

(1) ^ a + 2 2m{2mr - 2a - 3) 

2 ^ 2m2r2-2m(2a + 3)r + (a + 2)(2a + 3) ■ ^ ' 



The behavior of the potential (187) at A; = 4, m = 1 and a = 1 is shown in 
Figure 5. 
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Fig. 5. Scalar Coulomb potential So{r) (curve 1) 
and its Darboux transformed partner sl^\r) 
(curve 2). 
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10.3 Spherically symmetric potentials 



The first important detail to be taken into account is that the usual (3+l)-di- 
mensional Dirac equation with a spherically symmetric potential is included 
in our developments. Indeed, using the standard technique of separation of 
angular variables (see e.g. [40]), we get the following radial equation coupled 
to scalar Ws(r), pseudoscalar Wp{r) and vector W^ir) potentials: 



k 



(73 + [m + Ws\<T^ + i[E - W^]a2 ) ^{r) = 0, 



(189) 



where r is the radial variable, m is the mass of the particle, and E its energy, 
while k = ±1, ±2 ... is related to the total angular momentum. Eq. (189) can 
be also written as 




- + Wp 
r 



ai + [m + W,] a3-[E- W,] } ^{r) = , 



(190) 



which if Wy{r) = coincides with hoijj{r) = Eipi^r), ho = jdr + Vo(r), being 
Vo(r) given in the canonical form (19) with 



Po{r) = m + Ws{r) , qQ{r) = Wp{r) + 



k 



(191) 



Let us start with an unphysical potential, corresponding to po{r) = qo{r) = 
and A; = 0, i.e. Vq — mas. The associated Dirac Hamiltonian will be denoted 
by ho and its eigenfunctions by iIj{E), where we omit the evident dependence 
on the variable r and we stress only the dependence on the eigenvalue E. 
We will show that, starting with h^, one can obtain physically meaningful 
potentials with k ^ 0. For this purpose, we shall use the eigenfunctions of ho 
with E — ±m: 



I 1 I ~ ( —2mr 

iP[m) =1^1' = I ^ 



(192) 



^{-m) =1 , ^{-m) = I 

1 / I -2mr 



(193) 
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Other eigenfunctions will be also used for producing new nontrivial potentials 

^ cosh kr ^ 

(194) 



sinhA;r 



-kr 



I 1 \ 



\ A+m / 



m 



A2>0. 



(195) 



Example 12 Let us take the spinors Ui and U2 from which the transformation 
function lA is constructed as follows: Ui = ip{m), U2 = ipl—m). After a very 
simple algebra, we obtain the potential 

Vi = m(T3 + -(Ti, 

r 

which is just the free particle Dirac Hamiltonian with /c = 1. It is not difficult 
to see that the choice ui = ip{—'m), U2 = ip{m) gives the same Hamiltonian, 
but with k — —1. 

Now, all solutions of the Dirac equation with the potential Vi = m(T3 + -(Ji can 
be found either by applying the transformation operator (13) to solutions of 
the same equation with Vq = ma^, or by applying formulas (32) and (28)-(30). 
Two solutions with E — m are 



ip{m) 



1 

2mr 



ip{m) 



(196) 



For E — —m one gets 







<^(-m) =1 I , if{-m} 
1/r 




(197) 



The solutions with E ^ ±m have the form 



(/?(A) = LV'(A) = k 



— sinh kr ^ 
, cosh kr + ^ sinh kr , 

\ A— m (\—m)r / 



(198) 



<^(A) = L^(A) = ke 



-kr 



k t 1_ 
A Ar 



(199) 
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If we take now the spinors Ui = ip{m) and U2 = ip{—ni) for constructing the 
transformation function U of the second transformation step, we find after 
some algebra the free Dirac Hamiltonian with k = 2. 

In the next example we illustrate how the use of a spinor having one of its 
components equal to zero inside the transformation function produces a pseu- 
doscalar potential. 

Example 13 The spinors Ui — (p{m) and U2 — (p{X) give the simplest pseu- 
doscalar potential with the mass equal to —A 

?2(r) = — --4^, (200) 
r 1 + kr 



whereas the choice U2 = V^(A) corresponds to 

, , 2 3A;r cosh /cr — (3 + /c^r^) sinh/cr 



Another possibility to get a pseudoscalar potential is by choosing Ui = 'f{X) 
and U2 = By this means we get the potential corresponding to the 

same mass m and 

, . 2 k cosh kr + fc^r^cosech kr — 2 sinh kr , . 

*M = ; .(;^. cosh sinh ;^.) ■ (^°^' 



In the final example we are going to consider in this paper, a superposition of 
scalar and pseudoscalar potentials is shown. 

Example 14 We use here a linear combination of two spinors with eigenvalue 
m: u\ = apirn) + (p{m) (c is a constant) and U2 — (p{—m). Then, we obtain a 
potential with 



, , 16 (cm '^r^ + 2m^r'^) 

P2{r) = 2m- ^ ^ 3^3 , 203 

1 8m^r(2c + 3mr) 
q2{r) = -- + 8^3^3 + c(3 + Am^r'') ' ^ 

This potential is regular provided c > 0. 

Finally, we choose a transformation function composed by following spinors: 
Ui = ip{m) and U2 = V^(A), and we produce another exactly solvable potential 
corresponding to 
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P2{r) 



2kmr\ cosh kr + (A;^ + m(m — A)) sinh kr 



(205) 



2kmr cosh kr -\- {\ — m) sinh /cr 



q2{r)^k 



2kmr sinh /cr — A;(m + A) cosh kr 



(206) 



2kmr cosh /cr + (A — m) sinh kr 



The Darboux transformation for a generahzed Coulomb interaction is consid- 
ered in [41]. 



11 Conclusion and outlook 

In this paper only a small part of the properties of differential intertwin- 
ing (Darboux) operators for the one-dimensional Dirac equation has been 
considered. Wc have shown that some properties known for the case of the 
Schrodinger equation, like factorization of Hamiltonians, one-to-one correspon- 
dence between spaces of solutions of equations related with an intertwiner, 
supersymmetry, and determinant formulas for chains of transformations, also 
have their counterparts for the case of the Dirac equation. We have shown 
with numerous examples that this technique is as easily as in the case of 
the Schrodinger equation. Therefore, we hope that this paper will stimulate 
further investigations in this field. 

To put an end to this paper, we would like to enumerate some subjects that, 
from our point of view, are worth to be investigated in the future: 

(1) To get purely scalar spherically symmetric potentials. 

(2) To find conditions for a chain of transformations to produce regular po- 
tentials. 

(3) To find interrelation between differential and integral transformation op- 
erators. 

(4) To find transformation of scattering data such as transmission and reflec- 
tion coefficients, phase shifts. 

(5) To consider transformations of periodical potentials. 

(6) To apply this technique for describing the scattering of particles at high 
energies. 
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